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Stellar rotation: Evidence for a large horizontal turbulence and 

its effects on evolution 
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Abstract. We derive a new expression for the coefficient Dh of diffusion by horizontal turbulence in rotating 
stars. This new estimate can be up to two orders of magnitude larger than given by a previous expression. As 
a consequence the differential rotation on an equipotential is found to be very small, which reinforces Zahn's 
hypothesis of shellular rotation. The role of the so-called /^-currents, as well as the driving of circulation, are 
reduced by the large horizontal turbulence. Stellar evolutionary models for a 20 M© star are calculated with the 
new coefficient. The new and large Dh tends to limit the size of the convective core and at the same time it largely 
favours the diffusion of helium and nitrogen to the surface of rotating OB stars, a feature supported by recent 
observations. 
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1. Introduction 

Differential rotation usually creates turbulent motions due 
to shear in rotating stars. In a normally stable radia- 
tive layer, the turbulence is usually much stronger in the 
horizontal direction perpendicular to gravity than in the 
vertical direction (Zahn HT)92|) . The reason is that in the 
vertical direction the stable temperature gradient needs 
stronger forces to be overcome than in the horizontal di- 
rection where, in principle, no forces are opposed to the 
motions. An argument in favour of these intense horizon- 
tal meteorological-like turbulent motions is given by the 
study of the solar tachocline by Spiegel & Zahn (1992). 
The tachocline is the transition zone between the rigid ro- 
tation in the radiative interior and the external convective 
zone, where rotation varies with latitude. Spiegel & Zahn 
show that if the horizontal turbulence is intense, then the 
tachocline is very thin as supported by helioseismological 
observations. 

The global result of the transport of the fluid elements 
by the horizontal turbulence is represented by a coefficient 
of viscosity v^. This coefficient is a very important param- 
eter in the physics of rotating stars in several respects: 

-1. If strong enough, the horizontal coupling expressed 
by the coefficient v\ Y makes the angular velocity Q nearly 
constant on isobaric surfaces (cf. Zahn|1992). In this case, 
the angular velocity f2 is constant on shells and the ro- 
tation law is said "shellular" by Zahn. If this is the case, 
the equations of stellar structure are greatly simplified, 
because they depend on one coordinate only (which is not 



just the lagrangian coordinate M r , but which has to be de- 
fined in an appropriate way) . This enables us to keep a 1- 
D equations scheme for stellar structure (cf. Kippenhahn 
& Thomas IT575I Endal & Sofia 11575) 1 . In the case of dif- 
ferential shellular rotation, Meynet & Maeder ( 2000 ) have 
shown that the scheme usually employed is incorrect, but 
that a consistent 1-D scheme may still be defined. 

-2. The various mixing processes of chemical elements 
play a major role in massive star evolution (cf. Heger et al. 
125551 Heger & Langer l2000l Meynet & Maeder 12550)1 . The 
horizontal turbulence reduces very much the efficiency of 
vertical transport of elements by meridional circulation 
(Chaboyer & Zahn lTQ52|l . This enables us to understand 
in a consistent way why the vertical transport of chemi- 
cal elements by the circulation is much smaller than the 
vertical transport of angular momentum by circulation. 
This is a clear constraint which results from solar obser- 
vations (Chaboyer et al.m)95a, 19953), as wei l as from the 
observations of massive stars (Maeder & Meynet 2000 ) . 

-3. The horizontal turbulence was generally ignored in 
the treatment of meridional circulation or of shear mixing. 
However, recent developments (Chaboyer & Zahn [1992 
Maeder & Zahn IT5551 Maeder & Meynet I255T1 Briiggen 
& Hillebrandt 2001 ) show that the horizontal diffusion by 
turbulence may also intervene in the expressions of the 
transport of chemical elements by meridional circulation, 
of the circulation velocity, of the diffusion coefficient by 
shear mixing, of the heat transport, etc... Interestingly 
enough, the numerical convergence of the 4th order scheme 
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of differential equations expressing the transport of angu- 
lar momentum and meridional circulation appears to be 
sensitive to the value of diffusion coefficient of horizontal 
turbulence. 

We do not consider here the effect of the magnetic 
field (cf. Spruit 2002), which may also play a role in the 
transport of angular momentum. In Sect. 2, we examine 
the reasons which demand a new estimate of v^. In Sect. 
3, we derive a new expression for v\y and some numerical 
estimates. Sect. 4 provides a discussion of the results. 

2. Reasons for a new estimate of the horizontal 
turbulence 

The usual expression for the coefficient of viscosity due 
to horizontal turbulence and for the coeffcient D n of hor- 
izontal diffusion, which is of the same order, is, according 
to Zahn (1992, Eq.(2.29)), 



suggests that 



D h ~ v h = —r \2V0r) - aU(r)\ 

Ch 



(1) 



where r is the appropriately defined eulerian coordinate of 
the isobar (Meynet & Maeder ISOOU)) . Apart from the case 
of extreme rotational velocities, the parameter r is close 
to the average radius of an isobar, which is the radius at 
P 2 (cos , d) = 0, namely for i9 = 54.7 degrees. U(r) is the 
vertical component of the velocity of meridional circula- 



tion, V(r) the horizontal component, a 



1 d In r J Sl 



and Ch 



2 dlnr 

is a constant of order of unity or smaller. This equation 
was derived assuming that the differential rotation (as de- 
fined by the ratio in Eq.(4) below) on an isobaric 

surface be small compared to unity. Indeed, there are sev- 
eral difficulties suggesting us to reconsider the expression 
for v^: 

— The first reason why the above expression is not sat- 
isfactory has been given by Zahn (1992) and it is related 
to the way Eq.(l) has been obtained. If we write the dif- 
ferential rotation at a colatitude ■& as 



Q(r, i?) = tt(r) + Q(r, i?) with 
fi(r,i?) = fi 2 (r)P 2 (cos??), 



(3) 



where P 2 (cos z?) is the Legendre polynomial of second or- 
der, we find that the differential rotation is a constant (cf. 
Sect. 2.6 in Zaun lE^ . 



fi(r) 



Ch 

5 



(4) 



This ratio is obtained when we use a coefficient z^ n given 
by Eq.(l) together with the expressions for the horizontal 
transport of angular momentum. This ratio is smaller than 
unity, but, as noted by Zahn l|1992[) . there is no reason for 
the amount of differential rotation being constant with 
r. On the contrary, the importance of differential rotation 
should depend on the value of f n , because the stronger the 
horizontal turbulence, the stronger is the homogeneisation 
of the angular velocity on an equipotential surface. This 



Sh(r) 
Q(r) 



should decrease for larger tv Likely, 



we could also expect that the importance of differential 
rotation varies with the rotation velocity, since the hori- 
zontal turbulence is itself generated by rotation. 

- Another point is related to the numerical mod- 
els (Meynet & Maeder, EUDI Maeder & Meynet, EMU)- 
During the course of the evolution, some models indicate 
that the coefficient Z? n of horizontal turbulence, as given 
by Eq.(l), is not so much larger than the coefficient of 
vertical diffusion by shear, especially at low metallicity Z 
where U(r) is small (see Fig. 6 in Meynet & Maeder 2000 
and Fig. 2 in Maeder & Meynet I200I[) . This is not very 
satisfactory for the validity of the assumption of shellu- 
lar rotation. We may remark in this context that this as- 
sumption would be much better if the diffusion coefficient 
of horizontal turbulence would be larger. 

- We may also note that physically the horizontal tur- 
bulence results from the differential rotation, while the 
meridional circulation with components U(r) and V(r) re- 
sults from the disruption of the thermal equilibrium on an 
equipotential. These two phenomena are generally, but not 
necessarily, related. An example is the case of a uniformly 
rotating star. There we have no differential rotation, but 
a breakdown of thermal equilibrium occurs unavoidably. 

3. The dissipation and feeding of turbulent energy 

Let us firstly examine the rate of dissipation of the tur- 
bulent energy. As shown by Zahn l|I992J) . we may write 
the rate of viscous dissipation of the energy present in the 
differential zonal motions on an isobar as 

Se t (r,#)=iA l Lm#^StfJ , (5) 

per mass and time and for an interval of latitude 5$. 
Taking into account Eq.(3), we obtain 

5i t (r,#) = v h sm 2 # n 2 2 (r) ( ^f" 8 tl ) 5d 2 , (6) 



( 2 ) Set (r, 0) = 9 (r) sin 4 d cos 2 M-d 1 



(7) 



The rate of energy dissipation is proportional to the square 
of the amplitude f^M on an equipotential. It is zero at 
the pole and equator and maximum at /^(cos??) = 0. 

There is an excess of energy on an isobar due to the 
differential rotation described by Eq.(2) compared to an 
average rotation. The velocity of rotation v(r, $) on the 
equipotential of average distance r is given by 



v(r, "&) = r sin + r sin #fi2 (r) P 2 (cos ■&) 



(8) 



For an interval of latitude the difference of rotational 
velocity 5v(r, i?) due to the latitudinal differential rotation 
on the equipotential is 



5v(r,0) = r sin n 2 (r) dP ^ C0S ^ 5$ 

dv 

= -3r sin 2 $ cos $ fl 2 (r)6d. 



(9) 
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We express here only in the velocity difference due to 
the shear on the equipotcntial. The excess of energy 
<5-Ediff(y, i?) over an interval 5$ due to the differential ro- 
tation in latitude is 



^r 2 sin 4 ?? cos 2 fl|(r) <W 2 . 



(10) 



Now, this small excess of energy over an interval Sfi will 
be smeared out in a dynamical timescale St^m ■ 

Let us estimate this characteristic timescale. On the 
isobar, the differential rotation due to produces a shift 
Sip in longitude for two fluid elements located at a differ- 
ence 8$ in colatitude in a time interval St 



Sip = Sfl St 



so st 

av 



(11) 



The meridional circulation has an horizontal velocity com- 
ponent V, which is pointing toward the pole in the exter- 
nal layers where U(r) < 0. This is due to the Gratton- 
Opik term, which is the term — 2 ^Qp fQpik ITT>51|) which 
appears in the equation for U(r) given for example in 
Eq.(4.29) by Maeder & Zahn 1(159*5) 1 . Due to the aver- 
age density p at the denominator, the Gratton-Opik term 
is largely negative near the surface, thus it acts so as to 
change the sign of the circulation U (r) , making it rising in 
the equatorial plane and descending along the polar axis. 
In the deeper layers, one has generally U(r) > 0, which 
means that the circulation rises along the polar axis and 
descends in the equatorial plane. Thus, in this case V is 
pointing toward the equator. The shift in latitude is given 

by 



r S$ = V-^jr St , 
av 

and this leads to 



Sip = 



n 2 v f dp 2 

r \d<d 



[Stf 



(12) 



(13) 



The complex motion in t? and <p due to the differential 
rotation on the isobar will tend to smear out the lati- 
tudinal energy differences as discussed above. As a typical 
dynamical timescale, we take the time necessary for this 
differential motion in ip to perform n axial rotations. We 
may consider an average of Sip("d) over the star 



Sip 



n 2 v 



(tdis) 



dP 2 



sin *& d'd = 2mr 



Thus, we get for the characteristic timescale 



*diff 



5mr 



3 n 2 v 



1/2 



(14) 



(15) 



The numerical factor is of course rather arbitrary. The 
ratio of the energy excess (Eq.JTU|)) and of the rate of 



viscous dissipation (Eq. Q) is of the order of this timescale 
and we write 



SE diS (r, fl) 
Se t {r, tf) 



tdifi 



(16) 



Using Eqs.Q and (|1C)|> . we obtain for the coefficient of 
horizontal turbulence 



( 3 



V80mr 



r 3 n 2 v 



(17) 



We may also estimate i\ by dividing the square of a typi- 
cal lengthscale (of the order of r) by the diffusion timescale 
given by Ea. (|15|l . We obtain exactly the same functional 
dependence in (r 3 fl 2 V)^ , with a numerical coefficient de- 
pending on the chosen lengthscale. Studying conservation 
of the angular momentum by taking into account the hor- 
izontal variations fl 2 of rotation leads to the following re- 
lation (Zahn HMR Eq.(2.27)), that relates fl 2 and v h 



v h n 2 (r) = - Q(r)r[2V-aU] 
5 



(18) 



where a is the same as in Eq.(l). This is the expression 
discussed in Sect. 2 , which implies that, if is given 
by an equation like Eq.(l), the ratio =^y is a constant. 

Now, we eliminate fl 2 between the two equations (|17|) and 
JTSJ- This gives for the coefficient of viscosity due to the 
horizontal turbulence 



v h = Ar (rft(r) V [2V-a*7]) ! 

3 



with A 



400mr 



(19) 



For n=l, 3 or 5 A w 0.134,0.0927,0.0782 respectively. 
This expression can be written in the usual form = 
| I ■ v for a viscosity, where the appropriate velocity v is 
a geometric mean of 3 relevant velocities: 

— A velocity (2V - aU) as in Eq.(l) by Zahn fiMty . 

— The horizontal component V of the meridional circu- 
lation. 

— The average local rotational velocity rfi(r). This rota- 
tional velocity is usually much larger than either U (r) 
or V(r), typically by 6 to 8 orders of a magnitude in 
an upper Main Sequence star rotating with the average 
velocity. 



4. Numerical results 

4.1. Comparisons and orders of magnitude 

Let us compare the present value of to that given by 
Zahn l|1992J) . We get the following ratio from Eq.(l) and 

CB9, 



Vh (present) 
v\XZahn) 



Ac h 



( rQV 
\(2V-aUy 



(20) 
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The quantities V and U have the same order of magnitude 
The numerical models below show that typically V 
and (2V — all) « V, thus we have the following order of 
magnitude, 



v^iypresent) 



Ac h ( — 



(21) 



v\XZahn) 

Thus, we see that the ratio of the two estimates of the 
diffusion coefficient is equal to the power | of the ratio 
of the local rotational velocity to the horizontal velocity 
of meridional circulation at the considered level. Let us 
consider a 20 M Q star with an average rotation velocity of 
220 km/s at the surface. At the middle of the MS phase, 
the vertical component of the meridional circulation lies 
between 3 ■ 10 -4 and 3 ■ 10~ 3 m/s as shown by the models 
below (see also Meynet & Maeder l2000fl . Thus, we typi- 



cally have 



Uh(present) 



of the order of 10 2 (cf. Fig. 1). Thus, 



z^h {Zahn) 

our estimate of the diffusion coefficient of the horizontal 
turbulence is much larger than the coefficient proposed by 
Zahn lfTMSl Eq.(2.29)) as given by Eq.(l). 

Let us now estimate the degree of differential rotation 
corresponding to this value of v^. From Ea. H18|l . we have 
with Eq.JTnjl, 



n 2 

fi(r) 



1 ({2V 
5 A \ r 



aUf 



nv 



(22) 



There is of course no coefficient Ch in this ratio. 
Numerically, this is 1 /5 of the inverse of the ratio given 
by Ea. H20fl . in which the value of — 1 would be used. 
This results from Eq.(l) and Eq.(18) relating and J^- 
Thus, with the above estimates, we obtain a ratio of about 
= 2 • 10~ 3 . As the value of ia, obtained in this work 

O(r) 

is much larger than the value given by Eq.(l), we see that 
quite logically the degree of differential rotation on an iso- 
bar is much smaller. The present value of the coefficient 
reinforces Zahn's hypothesis of shellular rotation. We also 
notice that the ratio =2 2 - is larger for slowly rotating stars. 

This is quite a consistent feature, because is growing 
with the velocity of rotation. 



4.2. Test with the evolution of a 20 M Q model 

In order to examine the consequence of the new coefficient 
of horizontal diffusion, we calculate stellar models for a 20 
M Q with composition X = 0.705 and Z = 0.02 with the 
same physics as in our recent papers (Maeder & Meynet 
2001}. The initial rotation velocity is 300 km/s, which 



corresponds to average rotation during the MS phase of 
about 240 km/s. Several expressions and diffusion coeffi- 
cients will be discussed numerically below, let us briefly 
recall them. The vertical component U(r) of the velocity 
of meridional circulation velocity is given by 
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Fig. 1. Values of the various diffusion coefficient in the 
interior of a star model of 20 M Q at beginning of the MS 
phase, with an age of 8.577 • 10 4 yr. and central hydro- 
gen content of X c = 0.702. That is to say during the ini- 
tial non-stationary phase of convergence of the rotation, 
where the velocities U(r) are large. The new -Dh is that 
given by Eq.lO with A=0.079. The old D h is that given 
by Eq.(l). K is the radiative diffusivity. -D s hcar is the co- 
efficient of diffusion by shears. D c g expresses the diffusion 
of the chemical elements by meridional circulation with 
account of the effects of horizontal turbulence. When not 
specified, the quantities shown are those for the new D^. 



(23) 



P is the pressure, Cp the specific heat, Eq and are 
terms depending on the CI— and /i-distributions respec- 
tively, up to the third order derivatives and on various 
thermodynamic quantities (see details in Maeder & Zahn, 
1998}. The term Eq expresses the driving effects of merid- 
ional circulation, while the term E^ expresses the /i- 
currents which tend to inhibit the circulation. The term 

is very important numerically, its origin in this expres- 
sion is more complex than could be thought at first sight. 
This expression also prevents infinite velocities at the edge 
of semiconvective zones. The term Eq expresses the driv- 
ing of the circulation by the fluctuations of density due to 
the breakdown of radiative equilibrium. 

The diffusion by shear instabilities is expressed by a 
coefficient £> s hcar, namely 



U(r) 



P 



pgCpT [V ad - V + {<plS)V*\ 



D. 



shear 



A(K + D h ) 



! V M(l + 7 i) + ( V ad-V 1 .ad) 
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p 



I V dhir 



(V-v) 



(24) 



where / is a numerical factor equal to 0.8836, K is the 
thermal diffusivity and (V — V) expresses the difference 
between the internal nonadiabatic gradient and the local 
gradient (Maeder 2001)- There is also the coefficient D c s, 
which expresses the contribution of the meridional circu- 
lation and horizontal turbulence to the diffusion of the 
elements fZahn HM^ . 



D, 



cff 



\rU{r){< 
30-D h 



(25) 



while the transport of angular momentum by circulation 
has to be treated explicitely as an advection. More details 
on these various expressions, on the hypotheses leading to 
them and on their domain of validity can be found in the 
given references. 

Fig.f shows the diffusion coefficients at the very be- 
ginning of the MS phase. There, the situation is non- 
stationary during f-2 % of the MS lifetime, until the ro- 
tation has converged toward an equilibrium profile, (in 
reality a part of this convergence, but probably not the 
whole, may be achieved during the pre-MS phase). In this 
temporary stage, U (r) is usually much larger (about a few 
10 -2 m s _1 ) than later in the course of evolution, where it 
is only of the order of a few I0 -3 m s _1 (Meynet & Maeder 
2000} . We point out the much larger value of the new Dh 
with respect to the old one. With the new Dh, we see that 
-D G ff is rather small with respect to -D s h C ar, while with the 
old -Dh, the coefficient D c g would have been larger than 
-Dshcar everywhere, and in particular by several orders of 
a magnitude close to the core. As to D s h ea r, the effect is 
opposite, the new -Dh makes it bigger since K is replaced 
by K + Z?h, when the //-gradient is small. 

Fig. 2 shows the various diffusion coefficients near the 
middle of MS evolution. Interestingly enough, the star 
shows 3 cells of meridional circulation. At the interfaces 
located at =ffi = 0.535 and 0.950, the nulling of U(r) pro- 
duces a kink in the curves of -Dh, D s h ea r and D g. The 
outer cell is the Gratton-Opik cell, due to the lower den- 
sity in the outer layers. The main inner cell is the usually 
dominant cell where U(r) is positive. There the circulation 
rises along the polar axis and descends in the equatorial 
plane, (thus bringing angular momentum toward the inte- 
rior). The third cell close to the core is not a well under- 
stood one. It was already present in some curves of Fig. 4 
in Meynet & Maeder (f200f)|) . The velocities here are very 
small and slighty negative. We interpret this third cell as 
due to a change of the second derivative of the angular ve- 
locity fl, which influences the expression of U{r) as given 
by Maeder & Zahn (|1998|> . (We also remark a kick in the 
curve of D s h oa r at ^j- = 0.41; it is produced by variations 
of the nearly vertical gradient of /z in some regions.) 

Fig. 2 tells us a lot about the diffusion coefficients in 
the stellar interior and their effects: 

- As discussed above, the new Dh given by Eq.JTJJ is 




. v h L 



0.4 



0.6 



0.8 



M/M 



Fig. 2. Values of the various diffusion coefficients in the 
interior of a star model of 20 M Q at about the middle of 
the MS phase, with an age of 7.066 • 10 6 yr. Same remarks 
about the coefficients of diffusion as in Fig.f. 



larger by about 2 orders of a magnitude with respect to 
the old one given by Eq.(f ), the new value is not far from 
the thermal diffusivity K. 

- The new Dh brings some change to -D s hcar- In regions 
where the /x-gradient is negligible, the ratio J ^ 3h ° ar of 
the coefficients of shear diffusion calculated with the new 
and the old -Dh behaves like — ^P h . In view of the values 
in Fig. 2, this means that D s h ea r in the outer regions is 
increased only moderately, currently less than a factor of 
two. Comparisons with Fig. 6 by Meynet & Maeder (2000) 
confirms the comparable order of magnitude of .Dshcar- 

- When V M >> (V a d — V ra d), a situation which occurs 



close to the convective core, the ratio 
This means that -D s h, 



(new) 



behaves 



like D ^ new '. 



D BhBl ,r (Old) 

is increased by a fac- 
tor of 100 in the internal regions close to the core. Such 
a change should normally strongly favour mixing in the 
star, however this is not so much the case, because pre- 
cisely in the regions close to the core D s h ea r is very small 
due to the very steep fi gradient, which limits the shear 
diffusion as shown by Eq.J23J. Close to the core, D c s is 
generally similar or larger than -Dshcar (this was particu- 
larly the case when the low Dh given by Eq.(l) was used). 

- Contrarily to the case of -D sri ear, -Doff is reduced by an 
increase of -Dh, as is evident from Ea.l)25|l. This can also 
be seen from a comparison between the present Fig. 2 and 
Fig. 6 by Meynet & Maeder |20_00| . where much larger 
values of D e g can be seen. 

- Last but not least, the old -Dh was often of the same 
order as the old -D s h ca r in some parts of the star. This 
was not satisfactory, in view of the hypothesis of shellular 
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Fig. 3. Distribution of hydrogen in models during the MS 
phase when X c « 0.33 (cf. Fig. 2). The broken line shows 
the profile for a model without rotation. The continuous 
line shows the H-profile for a model with an initial rota- 
tion velocity of 300 km/s with new Dh, while the dotted 
line shows the H-profile for a rotating model with same 
rotation velocity and the old Dh. 

rotation as mentioned in Sect. 2. The new Dh, which is 
much larger than the new Z? s hcar (cf- Fig. 2), solves the 
problem and makes the hypothesis of shellular rotation a 
much better one as also indicated by Ea. H22l) . 

Thus we see that a larger horizontal turbulence makes 
-Dshcar larger and D c g smaller. The situation is complex, 
since the relative importance of these two coefficients is 
not the same throughout the star. D s hcar always domi- 
nates at some distance of the stellar core, while D e ff tends 
to dominate near the core, especially if Dh is small. In 
addition, the ratio of these two coefficients is changing 
during evolution, as seen from Fig. 1 and Fig. 2. Thus, a 
change of Dh affects the evolution of a star in a complex 
way. These new results now seem kind of very similar to 
Heger et al. (2000). In a rough summary, we may say that 
a larger Dh tends to reduce or contain the size of the core, 
since D e ff which is important close to the core is reduced; 
at the same time the spread of the processed elements (He 
and N) up to the surface is larger. 

A larger horizontal turbulence Dh also reduces the hor- 
izontal /Lt-gradients and thus it limits the importance of 
the so-called /x-currents introduced by Mestel H1953llT§6lfl 
see also Theado & Vauclair, 120011 and by Palacios et al, 
2002). Quantitatively, the term which expresses the 
/z-currents contains a term A = ^, as shown by Eq.(4.30) 

by Maeder & Zahn ifTMljl and A itself goes like (A,) -1 in 
a stationary situation (Eq.(4.40) in above reference). This 
establishes the relation of Dh with the /i-currents. 



In addition, the horizontal turbulence also contributes 
to smear out the temperature and density fluctuations on 
an equipotential and this reduces the effects driving merid- 
ional circulation. Quantitatively, this is expressed by the 
term containing Dh in the expression of Eq in Eq.(4.37) 
by Maeder & Zahn l|1998[> . Thus, globally a higher D h re- 
duces both the terms driving the meridional circulation as 
well as the term which inhibit this circulation. In the nu- 
merical example, we see that the values of U (r) mentioned 
above are generally smaller than those found by Meynet 
& Maeder <|2UUU}l . 

Fig. 3 illustrates the effects of the diffusion coefficients 
Dh on the internal distribution of hydrogen. The model 
with the new and higher values of Dh has a convective 
core and a surrounding H-profile which is close to that 
of the non-rotating model, in particular we see that the 
H-profile close to the core is much steeper than for the 
rotating model with the old Dh. In the outer layers, the 
H-content of the rotating model with the new Dh is lower 
than for the other two cases which means than mixing has 
been more efficient there. These properties arc quite con- 
sistent with our previous discussion. Indeeed, the higher 
Dh reduces the coefficient D g , which was the largest one 
close to the core. This prevents the growth of the core and 
creates the steep /t-gradient just above it. Now, the larger 
Dh makes D s h oa r larger outside the region of the very steep 
//-gradient and this favours mixing in the outer layers. As 
a consequence, the enrichments in helium and nitrogen 
at the stellar surface are higher. This explains the rather 
paradoxical result that the model with the higher Dh has 
a slightly smaller convective core and at the same time a 
larger enrichment in CNO processed elements at the stel- 
lar surface. 

Fig. 4 illustrates the tracks in the HR diagram. We see 
that the model with the new (and large) Dh has a turnoff 
inbetween that of the model without rotation and that of 
the rotating model with the old coefficient Dh given by 
Eq.(l). This is quite in agreement with the H-profiles and 
the values of the mass fractions of the convective cores, 
which are 6.8, 7.1 and 7.7 M© when X c = 0.33 for the 
model with zero rotation, for the model with rotation, 
with the new and the old Dh respectively. As well known, 
larger cores make MS tracks extend to higher luminosities. 
We notice however that the intermediate track with the 
new Dh is slightly bluer than an average of the other two 
tracks would suggest. This is due to the larger enrichments 
of the outer layers in helium, which reduces the opacity 
and makes the star slightly bluer and brighter. 

Fig. 5 completes this picture by showing the evolution 
with time of the ratio N / H of the nitrogen to hydrogen 
ratios at the stellar surface. We see that the surface enrich- 
ment in nitrogen of the model with the new Dh given by 
Eq. (|19|l is larger than the one obtained with old Dh given 
by Eq.(l). This is quite consistent with what we have just 
seen above in Fig. 3. The larger Dh makes D s h oa r larger 
in the outer layers and the transport of the new helium 
and nitrogen to the surface is more important. This obser- 
vational consequence is particularly interesting in view of 
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Fig. 4. The HR diagram for the MS phase of 3 models of 
20 M Q at Z = 0.02. The broken line (lower turnoff) is for 
a non-rotating model. The continuous line is for an initial 
velocity of 300 km/s with the new Z? n given by Ea. (|19(l . 
The track with a dotted line (higher turnoff) is for the 
same initial velocity with the old given by Eq.(l). 

the new results by Heap (j2~002|l , who has shown very high 
N/H enrichments in OB stars up to about 50. Future de- 
tailed comparisons considering carefully the mass, veloci- 
ties and abundances of OB stars are necessary to examine 
whether the models with the new -D sn oar are better sup- 
ported by the observations. 

5. Conclusions 

The following conclusions have been obtained here: 

- By expressing the balance between the energy dissi- 
pated by the horizontal turbulence and the excess of en- 
ergy present in the differential rotation on an equipoten- 
tial which can be dissipated in a dynamical time, we have 
found a new expression for the coefficient of diffusion _D n 
by the horizontal turbulence in rotating stars. This new 
coefficient is typically larger by a factor 10 2 than the one 
proposed by Zahn 1)1992[) . 

-The differential rotation on an equipotential is found 
much smaller so that the hypothesis of shellular rotation 
by Zahn (1992) is reinforced. 

-A higher horizontal turbulence reduces the importance 
of the /i-currents and also to a smaller extent the driving 
of the meridional circulation. 

-Numerical models show in agreement with a physical dis- 
cussion that due to the different effects of the horizontal 
turbulence on the shears and on the transport of chemi- 
cals by circulation, a larger Z\ tends to contain the size 
of the core and at the same time to favour the spread of 
the processed elements up to the stellar surface. 



Fig. 5. Evolution as a function of time of the abundance 
ratios Nj H at the stellar surface of the 3 models with 20 
M Q considered. The meaning of the lines is the same as 
in Fig. 3 and 4. The continuous (higher) line is that of the 
rotating model with the new D^. 

-The tracks in the HR diagram obtained with the new 
and larger D\ x for rotating stars are in agreement with the 
above effects. 

It will be interesting to further explore the conse- 
quences of the larger suggested here for other stellar 
masses and evolutionary stages. 
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